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ABSTRACT 

The bright uhrahiminous X-ray source (ULX), M33 X-8, has been observed several 
times by XMM-Newton, providing us with a rare opportunity to 'flux bin' the spec- 
tral data and search for changes in the average X-ray spectrum with flux level. The 
aggregated X-ray spectra appear unlike standard sub-Eddington accretion state spec- 
tra which, alongside the lack of discernible variability at any energy, argues strongly 
against conventional two-component, sub-Eddington models. Although the lack of 
variability could be consistent with disc-dominated spectra, sub-Eddington disc mod- 
els are not sufficiently broad to explain the observed spectra. Fits with a ~ Eddington 
accretion rate slim disc model are acceptable, but the fits show that the temperature 
decreases with flux, contrary to expectations, and this is accompanied by the appear- 
ance of a harder tail to the spectrum. Applying a suitable two-component model reveals 
that the disc becomes cooler and less advection dominated as the X-ray flux increases, 
and this is allied to the emergence of an optically-thick Comptonisation medium. We 
present a scenario in which this is explained by the onset of a radiatively-driven wind 
from the innermost regions of the accretion disc, as M33 X-8 exceeds the Eddington 
limit. Furthermore, we argue that the direct evolution of this spectrum with increas- 
ing luminosity (and hence radiation pressure) leads naturally to the two-component 
spectra seen in more luminous ULXs. 
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1 INTRODUCTION 

Ultra-luminous X-ray sources (ULXs) are point-like objects 
with high (>10'^® erg s~^) X-ray luminosities that are not as- 
sociated with an active galactic nucleus (AGN) or, indeed, 
the central regions of a host galaxy (see Miller & Colbert 
2004; Roberts 2007; Gladstone 2011). The nature of these 
objects has been the subject of much speculation, with CCD 
resolution X-ray spectroscopy playing a major role in ad- 
vancing our understanding. Although other missions have 
played an important part (e.g. ASCA detection of possible 
state transitions in ULXs, Kubota et al. 2001), these results 
have predominantly come from the XMM-Newton mission. 
Its first major advance was the detection of a soft excess in 
the spectra of many ULXs, with a temperature consistent 
with that expected for the inner edge of an accretion disc 
around an intermediate-mass black hole (IMBH; e.g. Miller 
et al. 2003; Miller, Fabian & Miller 2004). However, later 
studies showed that the second, harder component in these 
spectra turns over within the XMM-Newton bandpass, and 
so appears much cooler and optically thicker than the cor- 



responding Comptonisation media in Galactic black holes 
(Stobbart et al. 2006). This is inconsistent with the identi- 
fication of a sub-Eddington state for an IMBH, and more 
indicative of super-Eddington accretion onto a stellar-mass 
black hole (Gladstone, Done & Roberts 2009). The appar- 
ent divergence of the spectra of more luminous ULXs into 
two components (see Fig. 8 of Gladstone et al. 2009) can 
be interpreted in terms of the emergence of a radiatively- 
driven wind at super-Eddington accretion rates, with the 
outflowing material thermalising the underlying disc emis- 
sion to produce the soft spectral component as predicted by 
e.g. King (2004), Poutanen et al. (2007). The hard compo- 
nent is then produced within the photospheric radius, with 
its characteristic optically-thick Comptonisation signature 
either the result of a thick shroud of Comptonising elec- 
trons around the hot inner disc, or perhaps a change in the 
opacity of the outer layers of the hot inner accretion disc 
itself (Middleton et al. 2011). Such a model can explain the 
startling lack of variability seen in many of these sources 
(Hcil et al. 2009). In those few cases where large amounts of 



© 2009 RAS 



2 M. J. Middleton, A. D. Sutton & T. P. Roberts 



Table 1. XMM-Newton observations of M33 X-8. 



OBSID 


obs. date 


useful exposure 


ofE-axis angle 




net count rate 


total counts 


bin 






(ks, MOSl) 


(arcmin) 


(xlO'^-' erg em^^ 


(ct s-i, MOSl) 






0102640101 


2000-08-04 


8.2 


1.2 


1.80 1.84 


1.84 


30000 


Medium 


0102640301 


2000-08-07 


5.3 


13.9 


1.72 1.79 


0.69 


7500 


Medium 


0102641001 


2001-07-08 


8.7 


10.5 


1.63 -> 1.67 


0.88 


14500 


Low 


0102642001 


2001-08-15 


10.7 


14.1 


1.72 1.78 


0.61 


13000 


Medium 


0102642101 


2002-01-25 


12.2 


10.7 


1.75 -)• 1.79 


0.96 


22600 


Medium 


0102642301 


2002-01-27 


12.2 


8.7 


1.57 1.61 


1.05 


24300 


Low 


0141980501 


2003-01-22 


3.5 


1.1 


1.11 1.16 


1.11 


6800 


Low 


0141980601 


2003-01-23 


12.8 


14.1 


1.92 -> 1.97 


0.63 


16600 


High 


0141980401 


2003-01-24 


7.2 


12.7 


1.93 2.00 


0.77 


12200 


High 


0141980801 


2003-02-12 


12.6 


1.2 


1.41 -> 1.44 


1.33 


26800 


Low 


0141980101 


2003-07-11 


6.1 


10.6 


1.49 1.54 


0.88 


10600 


Low 


0141980301 


2003-07-25 















Notes: The observation date is given in year-month-day format. The observed flux, /x, the background-subtracted net count rate and 
the total source counts (sum of MOSl + MOS2) are in the 0.3 - 10 keV range. OBSID 0141980301 had no MOS observations and was 
excluded. 



variability can be seen, it is likely that the angle of observa- 
tion coincides with the edge of the turbulent photosphere, 
adding extrinsic variability to the X-ray signal (Middleton 
et al. 2011). 

Thus it appears that high quality X-ray spectroscopy 
of ULXs is consistent with many of the characteristics pre- 
dicted for super-Eddington accretion (e.g. Begelrnan et al. 
2006; Poutanen et al. 2007; Mineshige & Ohsuga 2011). 
However, there is a sub-class of ULXs with luminosities close 
to 10^^ erg s~^, that do not show the spectral inflection at 
~ 2 keV indicative of the two-component spectra in more 
luminous ULXs, including several objects in the sample of 
Gladstone et al. (2009, Fig. 8, top row). Their spectra are in- 
stead well-described by a single thermal component of emis- 
sion. Understanding the nature of this emission and how it 
fits into the broader picture of ULX X-ray spectra is impor- 
tant as it may provide a link between mass accretion rates at 
~ Eddington through to super- and hyper-Eddington rates. 
A prime method for investigating this is to observe the evolu- 
tion of X-ray spectra with luminosity for these objects. Such 
studies have been undertaken for a number of more luminous 
ULXs on the basis of data from Chandra (e.g. Roberts et 
al. 2006), Swift (e.g. Vierdayanti et al. 2010; Kong et al. 
2010) and XMM-Newton (e.g. Feng & Kaaret 2006; Feng 
& Kaaret 2009), yielding important insights including the 
demonstration that the temperature-luminosity relation for 
the soft excess is indicative of the presence of a wind-driven 
photosphere, rather than a standard accretion disc around 
an IMBH (Kajava & Poutanen 2009). However, spectral evo- 
lution studies for less luminous ULXs are less numerous, and 
can suffer from moderate data quality (e.g. four objects de- 
scribed by a simple multi-colour disc model in Kajava & 
Poutanen 2009). 

An obvious means of increasing our understanding of 
these objects is to study the brightest example, M33 X- 
8. This is the closest established ULX, long known as the 
brightest persistent X-ray source in the Local Group (Long 
et al. 1981; Gottwald, Pietsch & Hasinger 1987). It is also 
one of only two ULXs (with M82 X-1) whose flux regularly 
exceeds 10~^^ erg cm~^ s~^ in the XMM-Newton bandpass, 



and so has been the subject of many studies (e.g. Dubus, 

Charles & Long 2004; La Parola et al. 2003; Parmar et al. 
2001; Takano et al. 1994). It has been observed on several 
occasions with the EPIC instruments on XMM-Newton (see 
e.g. Foschini et al. 2004; 2006). However, due to the mod- 
erate statistical quality of individual observations it has not 
been possible to constrain changes to the X-ray spectrum 
between observations (Weng et al. 2009). Here we present 
a flux binned analysis that shows that the X-ray spectrum 
of M33 X-8 does alter subtly with luminosity, and has a lu- 
minosity dependence that is consistent with the picture of 
ULX spectral evolution in the super-Eddington regime. 



2 DATA ANALYSIS 

The EPIC instruments on-board XMM-Newton observed 
M33 X-8 on 12 occasions of varying durations during the 
early years of the mission (see Table 1 for details on each of 
these) ^. As the source is often not visible in the PN obser- 
vations (as it lies off the edge of the chip) we restricted our- 
selves to using only the MOS data. The data were initially 
cleaned using standard procedures. Good time interval files 
were created from the > 10 keV full-field light curves of each 
observation by removing periods with a raised background 
count rate, and used to filter subsequent data extractions. 

Where the source itself appeared in the outer MOS 
chips we used circular data extraction apertures of 45 arc- 
second radius for both source and background extractions, 
and in the three cases where the source is on-axis (OB- 
SIDs: 0102640101, 0141980501 and 0141980801 respectively) 
we instead use 30 arcsecond radius regions. We proceeded 
to filter the data using SAS VlO and extracted source and 
background spectra and lightcurves over standard event pat- 
terns (< 12) and flags (= 0). The background subtracted 
0.3 - 10 keV count rate is given in Table 1 and suggests 



New, deep observations were obtained in the summer of 2010 
that are not public at the time of writing. 
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0.001 0.01 0.001 0.01 0.001 0.01 
Frequency (Hz) 

Figure 1. Average, re-binned (by a geometrical factor of 1.2), 0.3 - 10 keV PDS with Itr error bars. The PDS are extracted from 
hghtcurvcs made from the combination of the finessed, background-subtracted, source lightcurves of each flux bin. In practice this has 
meant losing some of the shorter available segments in order to obtain a larger frequency bandpass. In the low flux binned data we 
use 401 bins/interval and 4 intervals/frame, in the medium bin, 525 bins/interval and 4 intervals/frame and in the high flux bin, 431 
bins/interval and 4 intervals/frame. The white noise level is shown as the horizontal dashed line and it is clear that there is no constrained 
variability power density in any bin. This is consistent with a model where the spectrum is dominated by stable disc emission. 



that there may be slight pile-up in OBSIDs 0141980501 
and 0141980801 (note that OBSID 0102640101 was taken in 
small window mode, and therefore does not suffer from pile- 
up at these count rates). However, inspection of the patterns 
using the SAS tool epatplot showed that this is not signifi- 
cant in the MOSl and only marginally significant for double 
patterns in the MOS2 observations and as such should not 
impact a combined spectrum where the small effect is con- 
siderably diluted. 

The crude individual spectra appear to be well de- 
scribed by a single component with a quasi-thermal shape. 
To obtain flux estimates for each observation, we model 
the MOSl and MOS2 spectral data with a simple thermal 
continuum (absorbed nthcomp) together with the CFLUX 
model in xspec, which provides an error range on the flux 
based on the data uncertainty. We achieve good or accept- 
able flts to each of the individual observation datasets and 
obtain the absorbed fluxes given in Table 1. From the ob- 
served distribution of fluxes we break the observations into 
3 classes: low, medium and high flux levels. Assuming a dis- 
tance of 817 kpc to M33, consistent with recent X-ray sur- 
veys of the galaxy (e.g. TuUmann et al. 2011), these bins 
correspond to observed luminosities of < 1.36 x 10'^^ erg s^^ 
(low fiux bin), (1.36 1.52) x 10^^ erg s"^ (medium fiux 
bin), and > 1.52 x 10^^ erg s~^ (high flux bin). Using 
the FTOOL ADDSPEC we proceed to co-add the respective 
datasets (spectra and response matrices) and obtain a MOSl 
and MOS2 dataset for each fiux bin. 

As our selection of fiux bin limits is somewhat arbitrary, 
we repeat the model dependent analyses reported in Section 
3 on rebinned data, where the highest fiux observation of 
the low and medium fiux bins was promoted into the next 
highest bin. In all cases we find behaviour consistent with 
the results reported below, indicating our results are not 
simply an artefact of the fiux bins used. 

The total aggregated counts for each observation are 
given in Table 1 and, whilst the dispersion in the co-added 
data will be affected by the weighted differences between 
the observations, the data quality overall is very good for 
the low and medium flux bins (totals of ~ 83000 and 53000 
counts respectively). The data quality for the high flux bin 



is somewhat poorer (--^ 29000 counts), but still a substantial 
improvement on most ULX datasets in the XMM-Newton 
archive, and sufficient to constrain the spectra and draw 
conclusions on spectral variability in M33 X-8. 



3 FLUX BINNED SPECTROSCOPY 

An initial inspection of the data using simple empirical mod- 
els demonstrated that M33 X-8 does indeed display spectral 
variability between the fiux binned datasets (which can be 
trivially confirmed by inspection of both Tables 2 & 3, and 
Figures 2 & 3). Here, we focus on how the spectra evolve 
with luminosity, and analyse this evolution in light of a range 
of physical assumptions and models. 

3.1 Sub-Eddington models 

Several authors have used the high luminosity of ULXs 
to claim the presence of an intermediate-mass black hole 
(>500 M0, e.g. Kaaret et al. 2003; Miller et al. 2003; Miller, 
Fabian & Miller 2004), which must imply sub-Eddington 
mass accretion rates. In black hole X-ray binaries such 
states are characterised by a two-component model of op- 
tically thick disc emission (Shakura & Sunyaev 1973) and 
Comptonisation (see the reviews of McClintock & Remil- 
lard 2006; Done, Gierlihski & Kubota 2007). We test this 
assertion by applying a model containing both of these 
components together with neutral absorption and a con- 
stant to account for the differences between the detectors 
(in XSPEC: constant*tbabs*(diskbb+comptt))^. Whilst 
Table 2 shows this model provides an acceptable fit to the 
data in all three cases, the characteristic temperature of the 
high energy tail is ~ 1.2 - 1.4 keV and so is inconsistent 
with the high coronal temperatures seen in sub-Eddington 

^ In these and all subsequent fits a constant component is utilised. 
It is fixed for MOSl, and the MOS2 value never deviates by more 
than ±5%. A Galactic column of 1.39 x 10^" cm"'^ in the direction 
of M33 (Kalberla et al. 2005) is assumed, and used as the lower 
bound in fitting the absorption component. 
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Table 2. Best fitting sub-Eddington models. 



Flux bin Low Medium High 



TBABS* (dISKBD+COMPTT) 

0.0691q;o24 0.080±0.006 0.081±0.011 

feTin (keV) 0.29l° °3 0.66±0.18 0.67±0.21 

fcreo^p(keV) 1.40l°:°« i-27l?,:l^ l-40lo:?9 

T 9.18to,lt > 10.14 > 8.33 

(d.o.f.) 689.8 (655) 654.9 (614) 521.9 (478) 

NuUP 0.17 0.12 0.08 

TBABS*KBRRBB 

Nh 0.052±0.003 0.072±0.003 0.067l„;oo5 

a > 0.999 > 0.971 > 0.991 

m 50+"-f'** Q 42+51.98 ^.2+5'iA5 

m '^■0^'-().29 ^-^^-0.39 "■^'^-0.36 

(d.o.f.) 1080.1 (655) 717.3 (614) 570.6 (478) 

NuUP 2.0x10-2-^ 2.0x10"^ 2.0x10-^ 



Notes: Best-fitting parameters for the two sub-Eddington mod- 
els used. The units of tfie foreground column, Nn, are lO'^'^ cm^'^. 
fcTin is the temperature of the inner edge of the accretion disc, and 
'i'^'comp and T are the temperature and optical depth of the Comp- 
tonising medium, a is the dimensionless spin parameter and m is 
the 'eflfective' mass accretion rate of the relativistically-smeared 
accretion disc (with the latter in units of IQl* g s~^), assuming 
zero torque at the inner boundary. The table also shows the values 
for and the number of degrees of freedom for each model, and 
the null hypothesis probability for this model being an acceptable 
fit to the data. 

XRB spectra (indicated by the unbroken power-law con- 
tinua out to > 100 keV, see e.g. McClintock & Remillard 

2006; although sec e.g. Zdziarski ct al. 2005 for the ULX- 
like behaviour of GRS 1915-1-105). We confirm this physical 
difference by fixing the plasma temperature at 50 keV in 
our models, and measuring the resulting change in fit qual- 
ity. In all three cases the fit is poorer with a hot corona, by 
Ax^ of 134, 42 and 15 respectively for one more degree of 
freedom. This shows that a cool corona provides a 99.98% 
improvement according to the F-test for the high flux bin, 
and a substantially more significant improvement for the 
other bins, compelling evidence that the corona does not 
appear similar to that in standard sub-Eddington states. 

Additionally, the Compton tail in a sub-Eddington state 
is highly variable on short timescales providing an unam- 
biguous test for such a model. We characterise the variabil- 
ity using the excess rms of the 3-10 keV lightcurve in each 
flux bin (binned to 250 s), i.e. the variability above the Pois- 
son (white) noise level of the lightcurve (see Edelson et al. 
2002) . We flnd that there is no constrained variability above 
3 keV with 3cr upper limits of <9%, <8% and <11% for the 
three flux bins respectively. This strongly argues against a 
two-component sub-Eddington model for the data. Indeed, 
we note that the combination of a hot disc (in 2/3 fits), and a 
cool, optically thick and invariant corona, demonstrates that 
this ULX cannot be described by the cool disc plus power- 
law continuum toy model previously used to infer IMBHs in 
ULXs. 

If the data were to be described solely by disc emis- 
sion then we would not expect there to be any variability 
on anything other than the longest timescales (Wilkinson 



& Uttley 2009). As we have substantial evidence against 
the two-component sub-Eddington model and are now in- 
terested in the variability behaviour of the emission as a 
whole, we extract the Fourier-frequency dependent power 
density spectra (PDS) over the full energy bandpass. Weng 
et al. (2009) report the lack of high frequency variability 
in the individual observations. Here however, we can ex- 
tend this to longer timescales by selecting the shortest seg- 
ment of continuous lightcurve and taking integer number of 
intervals of this length across the remaining observations. 
This can require certain finessing to obtain the maximum 
amount of available data and broadest frequency bandpass 
(as this goes from 1/ [longest available individual segment 
length] to l/[2x bin size]). In Figure 1 we present the av- 
erage power density spectrum (PDS - extracted using the 
FTOOL POWSPEC) for cach flux bin, normahsed to rms^ units 
following geometrical re-binning with white noise included. 
Quite clearly there is no constrained variability over any of 
the available frequency bins, consistent with the source hav- 
ing suppressed red noise as seen in a handful of ULXs (Heil 
et al. 2009), and the emission originating in the accretion 
disc. 

However, a thin disc is a poor description of the spec- 
tral data as the Wien tail is not broad enough to flt the 
high energy emission, leading to jjoor fits (x^ of 1441.5/658, 
931.4/617 and 650.7/481 for an absorbed multi-colour disc 
blackbody spectrum - diskbb in xspec - fit to each flux 
bin respectively). Disc emission that has been smeared 
by relativistic effects is considerably broader and so we 
also fitted the spectral data with the xspec model CON- 
STANT*TBABS*KERRBB. The data and best fitting models 
are shown in Figure 2, and the fit parameters are detailed 
in Table 2. Even in this extremely relativistically broadened 
case we still find the models poorly describe the data, with 
obvious excesses above the model at high energies. It would 
therefore appear that none of the sub-Eddington models can 
describe both the broad shape of the spectra and also the 
lack of high energy variability. 

3.2 Eddington to super-Eddington models 

Examples of Eddington accretion rates (i.e. accretion rates 
closely approximating those expected at the Eddington limit 
for an object) have been inferred in the outburst of several 
Galactic X-ray binaries, for example V404 Cyg (Zycki, Done 
& Smith 1999), V4641 Sgr (Revnivtsev et al. 2002) and the 
neutron star system Cir X-1 (Done & Gierlihski 2003). How- 
ever, ULXs are observed to be persistently luminous, which 
would suggest a closer analogy to a long-lived outburst sys- 
tem such as GRS 1915-1-105 (see Remillard & McCUntock 
2006), although this source is also known to undergo dra- 
matic state changes on short timescales (e.g. Belloni et al. 
2000; Middleton et al. 2006). 

Taking a theoretical approach, we would expect that, 
as the Eddington limit is neared, the properties of the flow 
through the disc would change as advection becomes more 
important as a method for removing energy from the flow 
(Mineshige et al. 2000; Abramowicz et al. 1988). We would 
also expect that, where the flow is highly illuminated, ma- 
terial would be driven off the disc surface by radiation pres- 
sure, with the massive outflow forming a photosphere above 
the disc (see Poutanen et al. 2007). This should modify the 
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Energy (keV) Energy (keV) Energy (keV) 

Figure 2. Flux binned X-ray data for the low, medium and high datasets (MOSl in black, MOS2 in red) folded with the best-fitting 
smeared disc model in blue (tbabs*kerrbb). The residuals to the best fit are shown below the best-fitting model and demonstrate that, 
with even a highly smeared disc, the broad shape of the data is not well matched. 



Table 3. Best fitting Eddington models. 



Flux bin 



Low 



Medium 



High 



kTin (keV) 
P 

(d.o.f.) 
Null P 



kTin (keV) 
P 



0.f72 



f.90 



TBABS'DISKPBB 

-t-0.011 nizis+0-Oii 
-0.004 '^•-^ ^"-0.010 
;L 42 +0.06 



",+0.13 
-0.08 

0.52±0.01 
699.10 (657) 
0.12 



0.05 

0.56±0.01 
667.5 (616) 
0.07 



0.022 
0.013 
+0.16 
0.07 
+0.01 
0.02 

532.1 (480) 
0.05 



0.151 
1.51 
0.56 



TBABS* (dISKPBB+COMPTT) 

0.160±0.010 
1.22±0.04 



n 11 =+0.023 



(keV) 



X2 (d.o.f.) 
Null P 



0.54±0.01 
< 141.76 
< 5.85 
691.0 (654) 
0.15 



0.93±0.02 
0.62±0.01 
1.671 
9.63 
651.4 (613) 
0.14 



0.68^ 



7+0.56 
-0.25 
+3.23 
-2.82 



0.099±0.03 

'^•"J-0.02 
a+0.04 
-0.01 
1 cr, + 0.33 
^•^^-0.19 

io.i9i^:J? 

521.6 (477) 
0.08 



Notes: Best-fitting parameters for the Eddington and/or super- 
Eddington models. Here p is the dimensionless index of the radial 
temperature dependence of an advection dominated slim disc (see 
text), and the other variables are as per Table 2. 



shape of the continuum emission via absorption, emission 
and scattering, producing a fully thermalised, blackbody- 
like component. Here however, the spectral shape of M33 
X-8 is unlikely to be described solely by a photosphere cov- 
ering a hot disc, as the peak temperatures of the observed 
spectra are implausibly high for a pure photosphere (c.f. 
King 2004). 



We therefore arrive at the model of an advection- 
dominated Eddington 'slim' disc (Mineshige et al. 2000). 
Several authors (Foschini et al. 2006; Weng et al. 2009) have 
shown this model to be a reasonable description of the mod- 
erate signal-to-noise data of individual observations, with a 
reported requirement for an additional soft power-law con- 
tinuum to high energies in some cases. Here however, we 
fit the flux binned data with an absorbed p-free disc model 
(C0NSTANT*TBABS*DlSKPBB in XSPEC), where the index of 
the temperature profile p is a free parameter (T oc 
for disc temperature T and radius R). Such a model has 
been used as evidence for slim disc spectra in ULXs, as a 
predicted change in the profile from p = 0.75 for standard 
discs to p = 0.5 for slim discs is indeed seen in some objects 
(e.g. Vierdayanti et al. 2006; although also see Gladstone 
et al. 2009). The best-fit parameters and their 90% confi- 
dence limits are shown in Table 3. The model provides a very 
acceptable description of the low and medium fiux binned 
datasets, and the peak temperature is constrained to drop 
with the increase in flux. The model is a poorer, although 
still marginally adequate description of the high flux data; 
however its parameters are not signiflcantly different from 
the medium flux bin. 

Although a slim disc by its very deflnition does not fol- 
low the same radial temperature profile as the standard thin 
disc (Shakura & Sunyaev 1973), we would not expect the 
temperature of the inner edge to drop with increasing flux 
unless the radius had moved substantially outwards (con- 
trary to what is seen in Galactic black hole X-ray binaries). 
Similarly, we would not expect the index of the radial tem- 
perature dependence p to increase signiflcantly with flux ei- 
ther. We can see from Table 3 that the temperature drops 
by more than 0.3 keV yet the p value increases signiflcantly. 



© 2009 RAS, MNRAS 000, 1-8 



6 M. J. Middleton, A. D. Sutton & T. P. Roberts 



T 1 — I I I I I 1 1 1 1 — 111 + 




^ — \ — I I I ii t i 




1 - - 

1 10 
Energy (keV) 



Figure 3. The ratio of the data in each flux bin to the best fitting 
slim disc model (tbabs*diskpbb) for the lowest flux bin. The 
panels show the residuals from: a) the low flux bin; b) the medium 
flux bin; and c) the high flux bin. The MOSl data is shown in 
black, and MOS2 in red. It is clear that, besides the expected 
increase in flux, the medium and high flux binned datasets differ 
in shape with an apparent deviation in the continuum above ~ 
2 keV. 



Thus the change here is contrary to what we would expect 
for an advection-dominated Eddington flow. However, the 
decrease in temperature of the peak of the emission could be 
explained by the inclusion of a flux-dependent cooling mech- 
anism. In particular, cooling could occur via the launching 
of a wind, and Comptonisation in the optically thick plasma 
of the outflowing material. In Eddington/super-Eddington 
flows, the launching of material is dominated by the radia- 
tion pressure (as the effective gravity goes as l-L/L^titi) and 
so one would expect to see its effect increase with luminosity. 
Hence, in terms of the spectral evolution, we might expect 
to see both the disc cool, and an optically-thick Comptoni- 
sation medium emerge, with increased luminosity. We illus- 
trate what the data shows in Figure 3, where we show the 
spectral residuals for each dataset compared to the best flt- 
ting model for the absorbed p-free disc model in the lowest 
flux bin (cf. also Table 3). This shows that, besides the ex- 
pected change in the flux level, we also see changes in the 
shape of the continuum with flux. In particular, a compo- 
nent with turnover in the 2-3 keV range appears in the 



medium flux bin spectrum, and appears more pronounced 
in the high flux bin. 

We attempt to model this by introducing a Compton 
tail to the model (using the COMPTT model, such that we 
employed a constant*tbabs*(diskpbb-|-comptt) model 
in XSPEC) with the seed photons drawn from the soft compo- 
nent in order to provide a physically limiting case and allow 
the spectral parameters to be constrained. In the case of the 
medium flux binned dataset the fit is notably improved by 
this additional component (Ax'^ of 16 for 3 d.o.f.s), and it 
is also marginally improved in the high flux binned dataset 
(Ax^> of 11 for 3 d.o.f.s). Although we would expect the 
slim disc emission to dominate at lower luminosities we also 
flt this model to the low flux binned data. We obtain only 
a slightly improved fit (Ax^> of 8 for 3 d.o.f.s) and, as ex- 
pected, the disc dominates the spectrum with only a small 
contribution to the higher energy emission made via Comp- 
tonisation. In order to realistically constrain the properties 
of the best-fitting two component model, we freeze each best- 
fitting component in turn and determine 90% error limits on 
the model parameters. The resulting fit parameters and er- 
rors are presented in Table 3. They suggest that, as the lumi- 
nosity of M33 X-8 increases, the soft component decreases in 
temperature and becomes less advection-dominated whilst 
the hard component becomes stronger and cooler. A plot of 
the evolution of the best-fitting models is shown in Figure 
4. As a consistency check we apply these best-fitting models 
to each of the individual observations of M33 X-8 in each 
bin and, in every case, find a good or acceptable fit to the 
data. 



4 DISCUSSION & CONCLUSION 

Assuming that the population of ULXs contain stellar mass 
black holes rather than IMBHs, their spectral evolution 
should provide the much sought-after explanation for the 
changing structure of the accretion fiow at and above the 
Eddington limit. This requires repeat observations of bright 
ULXs at different luminosities. Whilst the relatively large 
amounts of data from M33 X-8 have been analysed by other 
authors (e.g. Weng et al. 2009, La Parola et al. 2003, etc), it 
has not previously been flux binned in an overt attempt to 
identify spectral evolution. Using this method, we have been 
able to test models for the spectral behaviour of this source. 
We argue strongly against a sub-Eddington two-component 
description, particularly one describing ULXs as IMBHs, 
based on the shape of and lack of variability in the harder 
spectral component, that results from Comptonisation. The 
average PDS of each flux bin further constrained the source's 
properties and showed that, across the full energy bandpass, 
there is no constrained variability above the white noise at 
any frequency. Whilst this could imply that the spectrum 
is dominated by stable emission from a disc, we show that 
neither a standard nor a relativistically broadened thin (sub- 
Eddington) disc is a suitable description of the data. 

Given the source luminosity of 1-2x10^® erg s~^, Ed- 
dington mass accretion rates are easily satisfled by a black 
hole mass of ^lOM©, consistent with the masses of Galactic 
black holes. We flnd that an Eddington slim disc provides 
a good description of the data with constrained differences 
in the temperatures and radial emission profiles of the low 
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Figure 4. Best-fitting super-Eddington, two-component models to each of the flux bins. The blue line shows the total spectrum, with the 
green and red lines showing the disc and Comptonised spectral components respectively. Although the spectra are broad and can suffer 
from model degeneracies, the evolution of a slim disc spectrum alone from low to medium fluxes appears unphysical, and so this model 
presents the most likely physical description of the flux binned data. The changing properties of the disc and relative contributions from 
the two components is consistent with the increased flux leading to greater wind production, which in turn leads to increased cooling 
of the disc and hence an increasingly less advection-dominated solution. This is naturally accompanied by the emerging optically-thick 
coronal component as more material is driven out from the disc. 



and medium flux binned data. However, tliese differences are 
inconsistent with the expected characteristics of disc accre- 
tion, and instead imply the presence of an additional cooling 
mechanism. Subtle changes in the spectra of the medium and 
high flux binned states, compared to the low flux state, are 
well modelled by including a Comptonised component. 

M33 X-8, then, turns out to be a crucial object for fur- 
thering our understanding of accretion at Eddington and 
super-Eddington rates. At its lower fluxes it appears well 
described by a slim disc model alone; however as its flux in- 
creases, an optically-thick Comptonisation media becomes 
more important, in effect 'stretching' the disc-like spectrum 
by cooling the disc component, and simultaneously provid- 
ing a harder, upscattered, spectral component. This is cru- 
cial because it appears to reflect the initial stages of the di- 
vergence of the ULX spectra we see at higher luminosities, 
where the spectra are well-modelled by a cool, disc-like com- 
ponent and a harder, optically-thick Comptonisation com- 
ponent (Gladstone et al. 2009). If so, it is likely the result 
of the emergence of the expected radiatively-driven wind. 

We sketch a toy model of a possible physical scenario 
to explain both this initial wind emergence phase, and the 
higher luminosity photosphere phase in Figure 5. In the top 
panel, the accretion rate is only just entering the Eddington 
regime. Here, the soft component is the outer, advection- 
dominated slim disc, which extends down to some radius 
(given as R2 in Figure 5). Within this radius the accretion 
flow is highly illuminated and we see the emergence of out- 
flowing winds down to some radius close to the innermost 
stable circular orbit (ISCO) beyond which we expect the 
flow to become highly turbulent (possibly producing the ex- 
pected, as yet unobserved, high energy, optically thin com- 
ponent). The unbound plasma in the wind Comptonises the 
underlying hotter disc photons, producing the hard compo- 
nent in the spectrum. As the luminosity increases, winds 
may be driven from further out, whereas closer in, the wind 
may become increasingly mass loaded as bound material at 



lower luminosities may now be lifted from the 'surface' of the 
disc. This further cools the underlying accretion flow leading 
to a cooler peak in the hard spectral component. In addition, 
the soft component is now cut off at a larger radius and so 
peaks at a lower temperature as well as being less advection 
dominated (due to being further out in the gravitational po- 
tential of the BH). This scenario qualitatively matches the 
implied behaviour of M33 X-8. The key remaining issue is 
then how this scenario relates to the spectra seen for more 
luminous ULXs. We emphasise that in this toy model, as 
the X-ray luminosity increases we would expect to see the 
wind being driven from increasingly further out from the 
ISCO. Eventually this would leave little or no disc emission 
beyond the outer launching radius, as the very massive wind 
would extend over most of the hot (X-ray emitting) regions 
of the disc and fully thermalise the underlying disc photons, 
producing the soft excess seen in ULX spectra. However, we 
would also expect to see 'bare' disc emission from within the 
inner launching radius of the wind (the photospheric radius) 
where the local gravity is so low that the radiation pressure 
has 'blown' the excess material away, leaving an approximate 
hot thin disc. The thermal Comptonisation spectrum then 
originates either in a corona of hot electrons tightly bound 
to the upper layers of the accretion disc; or perhaps is the 
spectrum of the accretion disc itself, as its opacity would 
alter at the high temperatures within the photospheric ra- 
dius. This scenario is shown in the bottom panel of Fig. 5, 
and can explain the spectral properties of the more luminous 
ULXs seen to have clear two-component spectra (see Fig. 8 
Gladstone et al. 2009, also Middleton et al. 2011). In these 
cases the wind/photosphere emission is thermally decoupled 
from the inner disc emission and so the relative amounts of 
hard/soft emission are most likely degenerate in inclination 
angle and mass accretion rate. 

This model can also consistently explain the associated 
variability properties of ULXs as in all cases the spectral 
components in the X-ray bandpass are stable (although we 



© 2009 RAS, MNRAS 000, 1-8 



8 M. J. Middl 




Figure 5. Our toy model for ULX disc behaviour as the luminosity reaches (top panel) and far exceeds (bottom panel) its Eddington 
luminosity. The top panel shows how, at or near to Eddington, we expect the majority of the disc to be advection-dominated (slim disc: 
red), but as the flow is highly illuminated further in (shown as within the radius R2), we expect material to be lifted from the surface 
of the disc by radiation pressure (from the inner photospheric radius, Rl close to the inner stable circular orbit, ISCO) which cools the 
underlying flow. This outflowing, hot material Comptonises the hot disc photons, producing the hard component in the spectrum, whilst 
the soft component originates in the slim disc emission from beyond the outer photospheric radius (R2). As the mass accretion rate and 
luminosity increases we would expect winds to be driven from further out in the disc, and increased mass loading of the wind closer in. 
This produces a cooler hot component and a cooler disc component (as the inner edge of the 'naked' disc emission is now also further 
out), as is seen in M33 X-8. The bottom panel then shows what we think may be happening in those ULXs where the luminosity has 
increased further, to be substantially super-Eddington. Here, the radiation pressure is so intense that all the excess material within R2 
has been driven off leaving a bare, hot, thin disc. The outer launching radius now extends to cover all of the remaining X-ray emitting 
disc, producing the soft cool component in the spectrum in a thick photosphere, whilst the hard component is from the hot bare disc. 
We suggest that the varying amounts of each component evident in different ULX spectra (c.f. Gladstone et al. 2009) may be due to 
degeneracies in inclination and mass accretion rate between different ULXs. 



predict that at higher energies there is Ukely to be an op- 
tically thin, highly variable component as seen in Galactic 
black hole binaries). However, if our line-of-sight intercepts 
the launching region of the photosphere at high luminosities 
where the wind is in the outer rather than inner disc, then 
we expect this extrinsic variability to lead to the hard com- 
ponent being highly variable (e.g. NGC 5408 X-1: Middleton 
et al. 2011). 
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Hence, by examining how the X-ray spectral and tim- 
ing properties of the nearest ULX vary with source flux, 
we have shown that it is consistent with the picture of a 
source accreting at the threshold of the super-Eddington 
regime. As its luminosity increases, it begins to betray the 
signature of a possible outflow from its central regions. This 
may be the emergence of the outflow that appears to domi- 
nate the characteristics of ULXs at higher luminosities. This 
work shows the power of considering both spectral and tim- 
ing data, using the highest quality datasets available from 
XMM-Newton. Clearly, if we are to develop a deeper un- 
derstanding of how ULXs work, and further investigate this 
picture of objects dominated by a radiatively-driven wind, 
obtaining an increased number of similar datasets must be 
a priority for this and future missions. 
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